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Abstract 

The interquark static potential for heavy mesons described by a massive One Gluon 
Exchange interaction obtained from the propagator of the truncated Dyson- Schwinger 
equations does not reproduced the expected Cornell potential. I show that no formula- 
tion based on a finite propagator will lead to confinement of quenched QCD. I propose 
a mechanism based on a singular nonperturbative coupling constant which has the 
virtue of giving rise to a finite gluon propagator and (almost) linear confinement. The 
mechanism can be slightly modified to produce the screened potentials of unquenched 
QCD. 
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1 Introduction 



The resolution of Dyson-Schwinger equations leads to the freezing of the QCD running 
coupling (effective charge) in the infrared, which is best understood as a dynamical 
generation of a gluon mass function, giving rise to a momentum dependence which is 
free from infrared divergences [TJ[2] • Recently, we have calculated the interquark static 
potential for heavy mesons by assuming that it is given by a massive One Gluon Ex- 
change (OGE) interaction which we have called DS potential [3J. To our surprise the 
DS potential does not contain the physics of confinement in the quenched approxima- 
tion to Quantum Chromodynamics, namely the linear rise at large distances. In Fig. 
[H (right) the DS potential is shown for the parameters that fit the lattice propagator 
of ref. jl], using the mass and coupling constant equations of ref. [3] as I will recall in 
section 2, as seen in Fig. [T] (left). Note the finiteness of the gluon propagator. 

It is well known that a propagator with a functional form 1/g 4 leads to a linearly 
rising potential [5], however this propagator is singular at the origin contrary to the 
result of the lattice calculation of ref. [I] . Finite modifications of the confining Gribov 
propagator, ~ l/(<? 2 + m 2 ) 2 , have been proposed. These modifications lead to the 
description of chiral symmetry breaking through confinement with a mass paramenter 
~ 200 MeV Does this modified Gribov (mG) propagator describe confinement 
in the quenched approximation? It is clear that its funtional form in r space is a 
decreasing exponential and therefore in some range behaves linearly. Is the parameter 
space adequate to support linearity? 

I proceed by adding the mG propagator to the DS propagator [8] and study there- 
after the corresponding potential. I am able to fit the propagator data quite precisely 
with this Ansatz. However, the corresponding potential produces almost no confine- 
ment. The value obtained for the mass paramenter of the mG propagator is too large 
and the behavior is exponentail even for small values of r. If , on the contrary, I fit 
the Cornell potential, I get a Gribov type propagator, i.e. the mass parameter of the 
mG propagator tends to zero, leading to a huge rise close to the origin in disagreement 
with lattice QCD. 

The corollary of my mathematical analysis of the propagator data and the Cornell 
potential is that a good potential requires a singularity at the origin in momemtum 
space, while the propagator is finite in the physical region. One way out of the impasse 
is to introduce a singularity in the nonperturbative coupling constant. Redefining the 
potential in terms of this singular coupling (~ 1/g 4 ), I will analyze in detail two 
cases for which I am able to reproduce both the propagator and the potential: i) 
Dyson-Schwinger OGE and ii) mG + Dyson-Schwinger OGE. Finally, by softening the 
singularity I am able to fit a screened interquark heavy meson potential [3]. 

The results of this investigation will be presented as follows. In section 2, 1 introduce 
the formalism and discuss the changes after incorporating the modified Gribov term 
in the propagator. In section 3, I study the corresponding potential. In section 4, I 
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Figure 1: Left: Fit to the lattice propagator data of ref. [I] using the logarithmic 
mass equation and the following set of parameters for the Dyson- Schwinger propagator 
(solid) A = 0.300 GeV, 5 = 1/11, p = 2.0, m = 0.374 GeV, /i = 4.5 GeV, p x = 20.0, 
c = 0.27. Right: The "massive" One Gluon exchange potential for the same parameters 
(long-dashed) [3] compared with the Cornell potential (dotted). 



introduce, as a consequence of mathematical reasonings, a singularity in the coupling 
constant and discuss its consequences. Section 5 is dedicated to the formulation of the 
screened potentials and I finish in section 6 by drawing some conclusions. 



2 The Gluon Propagator 

Infrared finite solutions for the gluon propagator of quenched QCD are obtained from 
the gauge-invariant non-linear Schwinger-Dyson equation formulated in the Landau 
gauge of the background field method. These solutions may be fitted using a massive 
propagator. 

At the level of the Dyson-Schwinger equations (DSE) the generation of such a mass 
is associated with the existence of infrared finite solutions for the gluon propagator, 
A(g 2 ), i.e. solutions with A _1 (0) > 0. Such solutions may be fitted by "massive" 
euclidean propagator of the form 

A D 1 s (q 2 ) = m 2 (q 2 ) + q 2 (l + Tl(q 2 ,p 2 )) (1) 

where m 2 (q 2 ) and U[q 2 ,p 2 ) depend non-trivially on the momentum transfer q 2 . One 
physically motivated possibility, which we shall use in here, is the so called logarithmic 
mass running, which is defined by 

m 2 (q 2 ) = ml 
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Figure 2: Fit to the lattice propagator data of ref. [I] using the logarithmic mass 
equation and the following set of parameters: a) Dyson-Schwinger propagator (dashed) 
A = 0.300 GeV, 5 = 1/11, p = 1.0, m = 0.50 GeV, p = 4.5 GeV, p x = 1.0, c = 0.20 
; b) Confining propagator (dotted) s = 1.044 GeV 2 , m c = 0.390 GeV. The total 
propagator is represented by a continuous line. The inset reproduces the behavior 
close to the origin. 



where m , p and S are parameters whose values are chosen to fit the lattice propagator 
and A is the QCD scale. 

In order to fit the lattice data at a specific scale p 2 the following functional form 
has been proposed [5], 



1 + cln 



' q 2 + pi m 2 (q 2 ) 

To this propagator I add an effective propagator of the form, 



(3) 



sK 1 nf ( q 2 ) = (q 2 + m 2 c ) 2 , (4) 

where s is a dimensional constant and m c a mass determing the value of the full propa- 
gator at the origin. This behavior has been proposed to describe how confinement leads 
to chiral symmetry breaking through a gap equation [6j[7]. Thus our full propagator 
becomes 

A^(q 2 ) = <5^A(g 2 ), 
A(g 2 ) = A D5 (g 2 ) + A con/ (g 2 ). (5) 

In Fig. [21 I show how the lattice data can be fitted with this new propagator and the 
functional forms for the DS gluon propagator presented above 0. 

1 Note that in ref. [5] other functional forms have been used. 
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3 The Modified Dyson-Schwinger Potential 



Due to the presence of this dynamical gluon mass the strong effective charge extracted 
from these solutions freezes at a finite value, giving rise to an infrared fixed point for 
QCD. The non-perturbative generalization of a(q 2 ) the QCD running coupling, comes 
in the form 

' q 2 + pm 2 (q 2 Y 



a(q 



A) In 



A 2 



(6) 



where a = ^- and we take /3q = 11 — 2n//3 where n/ is the number of flavors. The 
m(q 2 ) in the argument of the logarithm tames the Landau pole, and a(q 2 ) freezes at a 
finite value in the IR, namely a -1 (0) = /3 ln(pm 2 (0)/A 2 ) [fl[2] . 
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Figure 3: The dashed curve shows the potential calculated using EqJH] with the pa- 
rameters obtained by fitting the lattice propagator of ref. [I] using the full propagator 
Eqs. E ■ A = 0.300 GeV, S = 1/11, p = 1.0, m = 0.50 GeV, s = 1.044 GeV 2 , 
m c = 0.390 GeV and particularizing for rif = 4, i.e. (3q = 25/3 in EqEJ The dotted 
curve represents the Cornell potential as given in ref. [3]. 



The potential between the heavy quarks is obtained from the One Gluon Exchange 
potential defined from the following propagator, [3] 

AnC F a(q 2 ) 
q 2 + m(q 2 ) ' 

where Cp is the Casimir eigenvalue of the fundamental representation of SU(3) [Cf = 
4/3]. To this term we add the modified Gribov propagator in Eq. HI 

The potential between static charges is related to the Fourier transform of the time- 
time component of the full gluon propagator which after some trivial algebra becomes. 
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V(?) = -A p|q| |q| (sC F - ^ + - * ) sin(|q||r|). (8) 
|r| Jo \ q z + m(q z ) (q^ + mcj^y 

For constant mass parameter the mG propagator has an exact Fourier transform leading 
to 

g— r m c 

where a = s/Ait. 

This potential for sufficiently small m c the potential behaves as linearly rising 
Kon/( r ) ~cr(— + | + . . .). Do the data support a small enough mass parameter 
to define a (almost) -liner behavior for the values of r required to fit the spectrum? 
In order to compare with the Cornell potential we have to implement the Sommer 
substraction as described in ref. [3j[9]. 

In Fig. [31 I compare the potential obtained from the fit to the lattice propagator 
of the previous section to the Cornell potential [T0HT3] . I note that the confining mass 
is too large to produce a linear rise at the relevant r values. Therefore the mechanism 
described above does not provide the required dynamics to describe quenched QCD 
confinement 1 1 1 . 

Let me now proceed in the opposite way. I fit the Cornell potential using Eq. [8] 
and construct the corresponding propagator. It is clear from the fit that the Cornell 
potential requires m c very close to zero. What happens then to the propagator? 




Figure 4: Left: Modified Dyson Schwinger potential (dashed) defined by the Slowing 
parameters: A = 0.300 GeV, 5 = 1/11, p = 2.0, m = 0.37 GeV, s = 4.316 GeV 2 , 
m c = 0.0 GeV and (3q = 25/3. The dotted curve represents the Cornell potential as 
given in ref. [3]. Right: The dotted curve represents the propagator obtained by using 
the parameters which fitted the Cornell potential. The continuous curve, which is 
almost indistinguishable from the data corresponds to the same fit without mG term, 
i.e. taking s = 0. 
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In Fig. HJ I show the fit to the Cornell potential and the corresponding propagator 
together with the lattice data. It is clear that the fitting of the Cornell potential 
requires a svery mall mass and a large coupling s, while the propagator requires a 
small coupling s if we choose a small mass. I show in the figure also the extreme case 
with s = 0, i.e. only with OGE propagator, which turns out to give a very precise 
fit to the propagator. Thus, one can fit the propagator with only the DS term as was 
done in ref. [8]. The corresponding parameters are those given in the figure caption. 

I am arriving to a impasse. The parameters which fit the propagator do not fit the 
potential and viceversa. The potential requires a confinement term a la Gribov while 
the propagator is finite at the origin. I therefore conclude from the analysis that it is 
impossible to fit the propagator and the potential simultaneously, since the potential 
requires a Gribov singularity for quenched QCD. 

4 Confinement Coupling Constant 

How can we solve this puzzle? The way I for see is dynamical. I need a singularity at 
the origin either in the confining mG or in the DS propagators in order to reproduce 
the Gribov behavior. However, the lattice data seem to imply that there is no singular 
behavior near the origin. One could still aim mathematically at a "hidden" singularity 
below the first data points. This solution does not make sense because a very narrow 
singularities in momentum space, becomes a fast flattening in coordinate space and 
therefore no linear rise will appear from such mechanism. My proposal here is to 
change the coupling constant and incorporate there the singularity. 
I have studied two mechanisms 

i) The Dyson- Schwinger scheme: the coupling constant changes as, 



This scheme has no mG term and the additional contribution could arise from a 
nonperturbative vertex correction. 

ii) The Gribov scheme: the potential acquires an additional nonperturbative cou- 
pling multiplying the modified Gribov propagator, 





A 4 



(10) 
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Figure 5: Left: Fit to the lattice propagator data of ref. [I] using the logarithmic mass 
equation and the following set of parameters: Dyson-Schwinger propagator (dashed) 
A = 0.300 GeV, 5 = 1/11, p = 2.0, m = 0.37 GeV, p = 4.5 GeV, p x = 20.0, 
c = 0.27. Right: The dashed curve shows the potential calculated using Eq{13]with the 
parameters obtained by fitting the lattice propagator, particularizing for rif — 4, i.e. 
0o = 25/3 in Eqj6] and using a con /(o) = 0.5. The dotted curve represents the Cornell 
potential as given in ref. [3]. 

The idea behind these Ansatze is that the strength of the 1/g 4 singularity, needed 
at the origin to achieve a linearly rising potential, cannot come from the propagator, 
which is finite, and therefore must come from the vertex. 

The potentials become now 



i) 

v(r) = -n fdlql l q l ( 8Cf JT' ( ?1U MlqlkD- (12) 
|r| Jo \ q + m (q ) J 

ii) 

V(r) = ~ p|q| |q| UC F " (q2) + y(q 2 ) ^ sin(|q||r|). (13) 
|r| Jo V + m(q z ) (q 2 + m c ) z J 
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Figure 6: Left: Fit to the lattice propagator data of ref. [I] using the logarithmic mass 
equation and the following set of parameters: a) Dyson-Schwinger propagator (dashed) 
A = 0.300 GeV, 5 = 1/11, p = 2.0, m = 0.37 GeV, fi = 4.5 GeV, p x = 1.0, c = 0.2; b) 
Confining propagator (dotted) s = 0.094 GeV 2 , m c = 0.390 GeV. Right: The dashed 
curve shows the potential calculated using EqJTjJ] with the parameters obtained by 
fitting the lattice propagator, particularizing for rif = 4, i.e. (3 = 25 /3 in Eqj6] and 
using a con /(o)S = 18.2 GeV 2 . The dotted curve represents the Cornell potential as given 
in ref. [3] . 

In Figs. |5] and El I show the result of the above procedure for the two cases. 

The fits in the two cases are very good. In particular, for the DS case, the fit is 
almost perfect. I conclude that, as far as the propagator is finite, one way to achieve the 
dynamics of quenched QCD confinement in the scheme is by a singularities appearing 
in the vertex. 

5 Screened Potentials 

It is well known that in real QCD the coupling constant is finite [15]. In my scheme 
this can be achieved avoiding the singularity by introducing a cutoff mass into the 
coupling constant, i.e. q 2 — > q 2 + m 2 . From the point of view of interquark dynamics 
it is known that once the theory is unquenched, the linearly rising potential flattens 
at large r leading to the so called Aachen potential [3j[T6l[TT] . By using this cutoff 
mechanism I show the resulting potentials in Fig. [7J Recall that the propagator fits 
are the same as before since the Gribov mechanism does not affect them. Again the 
fit is excellent in both cases, almost perfect in the DS case. 
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Figure 7: The dotted curve represents the Aachen potential as given in ref. [3]. Left: 
Modified Dyson Schwinger with soft confinement potential (dashed) defined by the 
Slowing parameters: A = 0.300 GeV, 5 = 1/11, p = 2.0, m = 0.37 GeV, m c = 0.390 
GeV, a con /( ) = 0.5,m a = 0.105 GeV and (3 = 25 /3. Right: Modified Gribov with 
soft confinement potential (dashed) defined by the Slowing parameters: A = 0.300 
GeV, 5 = 1/11, p = 2.0, m = 0.378 GeV, a conf(0) s = 18.2 GeV 2 , m c = 0.390 GeV, 
m a = 0.090 GeV and /3 = 25 /3. 



6 Concluding Remarks 

The gluon propagator in the Landau gauge in quenched lattice QCD is finite. This 
has been shown to be the case also for the resolution of truncated Dyson-Schwinger 
equations. This finiteness implies a limiting value for the corresponing interquark 
potential which does not correspond to the potential obtained from quenched QCD 
which is linearly rising. I simplest solution I forsee is to implemant a Gribov singular- 
ity in the nonperturbative coupling constant. In my calculations all the parameters, 
except those related to the singularity, have been fixed to the lattice propagator and 
the correct potential arises from adjusting the singularity. In this way I am able to 
reproduce the Cornell potential with great precision. The screened potential, corre- 
sponding to unquenched QCD, arises naturally by modifying the Gribov singularity 
with an additional mass parameter. At this point we are not able to adscribe physical 
meaning to the parameters since they arise not from fundamental equations but from 
parametrizations. 

The conclusion of this study is that one is able to reproduce the gluon propagator 
and the heavy interquark potential if one is able to associate an interesting dynamical 
content to the nonperturbative coupling constant. A more fundamental study of this 
element from the point of view of nonperturbative QCD studies will shed some light 
in the confinement mechanism. 
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